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ABSTRACT 


A  precision  formula  is  here  given  for  the  self  inductance  of  a  single-layer  helix 
wound  with  ordinary  round  wire.  This  is  neither  more  nor  less  accurate  than 
that  previously  published,  from  which  it  has  been  derived  by  evaluating  certain 
correction  terms  and  replacing  them  by  approximation  formulas  which  are  much 
simpler  to  compute. 
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I.  INTRODUCTION 

The  formula  derived  in  this  paper  is  a  modification,  without  loss  of 
precision,  of  the  expression  for  the  inductance  of  a  single-layer  helix 
of  round  wire,  given  in  Bureau  of  Standards  Scientific  Paper  No. 
537  (vol.  21,  p.  431,  1926-27).  The  notation  has  been  changed,  and 
the  correction  terms  transformed  so  that  the  present  formula  is  more 
simple  from  the  point  of  view  of  the  computer. 

The  diameter  of  the  wire  is  d;  the  mean  diameter  of  the  solenoid  is 
D;  and  its  length  I  is  the  axial  distance  from  the  center  of  the  wire 
at  the  begininng  of  the  first  turn  to  the  center  of  the  wire  at  the  end 
of  the  Nth  turn,  N  being  the  total  turns  so  that  the  pitch  of  the  winding 
is  l/N.  The  modulus  k  of  the  complete  ellipitic  integrals  K  and  E  is 
given  by 

T)2 

k2  —  — =—— 

~~l2  +  D2 
The  principal  term  Ls  is  the  current-sheet  formula  of  Lorenz. 

As  one  of  the  correction  terms  there  appears  M,  the  mutual  inductance 
between  the  two  end-circles  of  the  solenoid,  which  is  computed  by 

M^D[^^k]  (2) 

The  inductance  L  of  an  actual  helix  is  (to  a  precision  which  neglects 

/   I   V         I 
terms  of  the  order  of  (  j^  )  log  t^) 
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L=Ls+,D[2N(ln±-0.89myiln^yip(£\ 

~\M-^F+W{E-k)[2±(^f)-\  (3) 

^["i      I    ■  _,       D     1 

2L     #       -y/p+ipj 


where 


(4) 


P(v)  =  j~i  +  2ln4:t]  when  ^  1 

—  3.if  —  17/71*7  when  77  ^  1 

The  ambiguous  sign  is  plus  for  the  "natural"  and  minus  for  the  uni- 
form distribution  of  current  over  the  section  of  the  wires.  A  current 
density  inversely  proportional  to  the  distance  from  the  axis  of  the 
solenoid  is  called  the  natural  distribution 

The  formula  for  L  here  given  is  obtained  by  change  of  notation 
and  by  transformation  of  the  function  A2  (k)  in  formula  (114)  of  the 

paper  cited.     (This  formula  contains  a  misprint  in  the  term-i  W  " 

la  3       P 

which  should  be  g  log  -) .     In  the  present  notation  this  formula  is 

+  7rDA2(k)  t  I  ^F+D\E-k)  f^\  (A) 

The  function  A2(k)  is  defined  by  formula  (89) 


Mfc)=-2[B0(k)-B1(k)+B2]  +  0.66267 -|^ 

o 

The  numerical  quantity  B2  was  denned  in  (90)  as  a  series  which  has 
been  found  to  be  -0.60835.    Hence 

?r2M2(£)  =  Z>£j(£0(*)  -  £,  (£))  +  o.i  il  (5) 

The  term  0.11  in  the  parenthesis  will  be  omitted  since  0.1  D  will  be 
considered  negligible  in  this  small  correction  term.    As  mostlole 
noids  to  be  used  as  precision  standards  would  have  an  Sductance 

ssBsart3S=3saa?5Sae which  '*- ,= 

II.  TRANSFORMATION  OF  B0(k) 
On  page  507  of  the  paper  quoted,  BQ  is  given  by 


(7) 
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If  we  let 

^Vl-fc2  and  z2  =  —L-2  and     . 
k  1  +  X2 

j>t  v-±  r  K(&d*  _4  r  d*  tc(  l  \ 

-±  fide  fi-rfi^r^-  (hog[r^X^]de 

xjo       Jo  ^\2  +  sm2d    ttJo      SL         sind        J 
then  (6)  gives 

;  ^BQ(k)=-^^¥+DnE(k)-k]  +  lpQ)  (8) 

We  may  evaluate  P(tj)  first  for  the  case  where  77=!;  that  is,  where 
1>D.     Using  the  formula 

2  log  sin  0  d0  =  - ^  log  2 
we  may  write  (7)  in  the  form 

•jr 

P(tl)  =  2  log  2  +-  f2log  (77+  V?+sm^)  d0 
ttJo 


or  since  (for  77=^1) 

logV     V     ^>V      s      2^Tr(»+i)T 

00 

PM     „,       .     ,    1      V(-l)Tl-3.5  ■  ■   (2w-l)7l    .,    _, 
Pin)  =2  log  4»>  +  i?-2J-irL2.4.6  .  .  .     2n    J  ^  lf  "£1 


(9) 


Since  the  sum  of  this  alternating  series  is  numerically  less  than  the 

9  1 

value  of  its  first  term  tk^~v  which  is  never  greater  than  about  ti^V 

g£l,  it  is  sufficient  to  take 

PW  =  2  log  4^+^         when  r^l  (10) 

When  i7  =  l  this  becomes  2  log  4  +  2  =  3.02  while  the  exact  expression 

(9)  becomes  2.98.     We  may,  therefore,  take  P(l)  =  3.     On  the  other 
hand,  when  0  =  77 =  1,  the  approximation 

P  (77)  =377  +  77  log  -         whenO^ry^l  (10)' 

V 
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will  be  in  error  by  not  more  than  0.1.     This  may  be  seen  from  the 
0  and  1  is       1S  C°rre      at  V==°  and  V  =  1'  While  the  error  between 

Y(V)==Sr]  +  Vlog  1~P(V) 

which  has  its  greatest  numerical  value  when 


Vl-P 
*=?£ —  =  0.34904 

corresponding  to  ¥  =  0.8914.     This  is  the  root  of  the  equation 
rM-O-2  +  log  ^-P'«  =  2+l  log  i-f  ^ 

ll^t^t^iS^ error' we  may  obtain  an  upper  and  Wr 

Write 

P(v)  =  log  4  +  -J     log  (r,  +  Vl  +  i?2-cos20)2d0 
1  f  T 

=  l0g  4+-Jo     log  (1+27?2-  COS2* +  2^1+**- COS2*)  de 

Now 

2t7VT+?:=^oW=2^Vi+72--  jj2=  cos20-Z  (,,  0) 


where 


CO 

V?r   Zjr(s+i) 


where 


The  function  Z  is  never  negative.    It  vanishes  when  fi=|  and  has  its 
greatest  value  when  6 = o  and  0  =  x  which  is 

Hence  P(ij)  may  be  written  in  the  form 

P(v)=lf;  log  ^JTR+ny-Zir,,  ti±(£&2+*)  cos^ 

IL%  *fis  V1^*1  we  rePlace  the  variable  function  Z(„,  9)  bv  its 
SSB  ft  K  f  a  fUnCti°n  ™'  Whi0h  iS  -  ^SJW 
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i  If  we  replace  Z(i?,0)  by  its  greatest  value,  we  get  a  function /(ij)  which 
is  a  lower  limit  to  P(y)  for  all  values  of  (ij)  where 


•de 


/(,^j>4((VI^,,-'i^x^)c„.,l 

That  is/(0)  <P(v)  <F(v)  if  ij>0. 

Both  of  these  functions  may  be  evaluated  by  means  of  the  definite 
integral  formula 

M^log  4  (a2-b2cos29)  dd  =  2  log  (a+  ^a^b2)  if  a>b 
which  gives  

[  i^)  =  2log  (Vl+^2  +  ^)(l+^/^=-2) 


This  gives 


/(n)  =  2log^2n  +  ^l  +  4i|J,  +  -T=L^) 


F(0)  =  0,F(1)  =2.984 

/(0)  =  0,  /(1)  =  2.958 

At  the  point  of  maximum  error  in  the  range  0<?7<1,  which  is  at 
77!  =  0.349,  we  find 

^(r?1)  =  l  43llwhence  p0?i)  =  1-506  ±0.075 
and 

3  rh  +  Th  log  —  =  1.415    • 

which  shows  that  the  maximum  error  made  by  using  (10/  for  P(t])  is 
approximately  0.1  in  the  range  0=  77  =  1. 

III.  TRANSFORMATION  OF  Bx(k) 

The  function  Bx(k)  is  defined  in  the  paper  quoted,  on  page  453,  by 

kB,(k)  =  ^  E(k)-  jl  d2jl-k2sm29dd 

(h2-y/l-k2sm2ddd^E(k)-  fh (*■- 0) V 'l -k2  cos2  <j>d<}> 
and 

Hence 


Now 
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where 
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™*-ff! 


*2    _V!^F^?0^ 


The  factor  of  the  integrand  y(#)-|,0  (.-*)  is  represented  by  a 
parabola  which  vanishes  when  *  =  0  and  has  the  maximum  value 
unity  when  #--,  where  it  also  has  a  horizontal  slope.    It  is  therefore 

SSFSffi&T*  "  SiU  *    *******  *  *  Sin  *  i*S  the 

£/(&)  =  _ ^  ft.  cos20sin0 
2  Jo  Vl~ 


*^*-^t-^ 


-P 


/j 


sin-1& 


o 


Hence 


-1 


£> 


yz2+z>2 

b        2L         Z>  -y/l2  +  D2J 


.] 


-r^ 


Figuee  1.— Arrangements  of  lead 

^^^^^i^^fk  **  mutual 
equations  (10)  or  (lOV ITfrn a.?lrcIes  °*  the  solenoid.  By  use  of 
formula  (3)  of  this  paper      (  **  ^  6q™tl0n  (A)  B  ^formed  into 

IV.  CORRECTIONS  DUE  TO  LEAD-IN  WIRES 
The  inductance  of  the  helix  and  lead-in  wires  together  is 

L'=L+Ll  +  2M„ 

where  L  is  that  of  the  helix  given  bv  formula  li\    r    ■    *i.       « 
inductance  of  the  lead-in  wires  anrf^^hi    (  I'  Z,'    S  ,the  self- 
between  these  wires  and  thThelix     In^omnuC  /^  m^T 
may  be  considered  a  current  shee",  £  ^^Jfc*^£ 
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of  current  contributes  to  Mhi.  If  the  leads  are  disposed  as  shown  in 
Figure  1,  it  is  evident  that  these  leads  which  are  perpendicular  to  the 
axis  of  the  cylinder  have  no  mutual  inductance  with  the  cylinder. 
The  latter  carries  unit  current  in  the  y  =  direction  so  that  the  current 

density  is  ^ —  =  — pr     Consequently,  the  term  Mhi  is  given  by 

Cb/2  Cl/2 

Mhl  =  2j0A(r,  y)dy-2j0A(r,  y)dy  (12) 

A  (r,  y)  =  2i  \  dy'  I  ^(y-y')2  +  a2  +  r2-2arcos<f> 

J -1/2    Jo 

.  »w    r*12  d  e 


where 


=  ~  I  dy'  j    -y/iy-y'y+ia  +  ry-lar  sin2  d 


-1/2    i/o 


2  f  I  f 

t/0  t/0 


V^2  +  (a  +  r)2  —  4  a  r  sin2  0 


+  -  [  dz        V22+  (a  +  r)2-4  a  r  sin2  (9 


0  t/0 


or 


ay  fcj  t/  &-y 

where 

V+7/  (14) 

As  a  special  case  of  this 


4  I  K  (x)  d  x     -,->  /    \  i      /«\ 


h2  = 7-t^tj  = /   7  m  and  ij, 


(15) 


J_V     1+/  !  \—"     2D  (16) 


Where 


,4  a  J      x  '  \2D 
Now  since  6  and  r  —  a  =  h  are  both  small,  we  may  write 

fb/2  rb/2 

2    A(r,y)dy  =  2A(a,Q)     dy  =  bA(a,0)  =  bP(Vo)  (17) 

Jo  Jo 

In  the  first  integral  of  (12)  we  may  write 

A(ry)  =  A{a  +  h,y)=A(a,y)  +  h^-^^=a 
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Computing  p*M]      from  (13),  since  ^«     _0 
L     °r     Jr=a  L>a  +  rJr=a~U' 

one  obtains  after  placing  r  =  a 

or  u  or  D 

so  that 

where 

tC    y  = 


1  +  (lJ2  +  y 
Changing  the  variable  of  integration  from  y  to  k,  gives 

=  _±D  f1  _dz_  C1  K(x)dx         4L> 

*  Jk  zyT=^Jz  x~/T^=  -~TBo(k)   (19) 

Adding  (17),  (18),  and  (19)  gives  by  (12) 


where 


¥  =  ~HF  and  *»2=  — iy>™ 
1+i  1+(A) 


By  use  of  (8)  this  may  be  put  in  the  general  form 

M1+5)v^[*-*]-*p(4)+»p(^)     (21) 


where  P(z)  is  defined  by  (4) 
Washington,  June  17,  1932. 


